Correlating elasticity and crack formation 
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For solving the longstanding materials science problem of correlating elastic properties of a solid 
material to the formation of cracks we present a new general concept. This concept is applied to the 
technologically most important cracks of loading mode I for which we establish exact correlations 
by introducing a localization length as a new material parameter. We study two limiting cases of 
crack formation making use of analytic models determining the material and direction dependent 
parameters by comparison to ab initio density functional results. This is done for a variety of real 
materials in order to test our approach for different types of bonding. For the most interesting 
ideal brittle cleavage we find that the localization length is -within a reasonable approximation- of 
constant value, which results in a simple relation for the critical stress, presumably being useful for 
materials engineering. Our results confirm that the proposed general concept results in meaningful 
physical models which -for the first time- allow a rigorous and simple correlation between elastic 
and mechanical properties. 
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1. INTRODUCTION 

It is a long standing effort in materials science to de- 
scribe or just to estimate the critical property of crack 
formation of a solid material in terms of its elastic proper- 
ties, an objective both of scientific as well as technological 
interest. Although more than 80 years have been passed 
since the pioneering work of Griffith [lj, it seems that 
still there is a need for quantitative models on the atomic 
scale in order to establish a correlation between elastic 
and mechanical properties. First attempts for cracks of 
mode I estimating the critical cleavage stress were made 
by Polanyi and Orowan Further studies tried 
to improve Orowan's approach but the quantitative re- 
sults were not satisfying. Very recently, Hayes et al. Q 
-based on the work of Nguyen and Ortiz jg- applied a 
rather complicated concept, which for intrinsic proper- 
ties on the atomic scale involves a dependency on the 
macroscopic thickness of the material. 

The conceptual problem of correlating elastic and 
cleavage properties consists in correlating a non-local 
property to a local property: the elastic response to a 
perturbation is non-local because the energy of is dis- 
tributed over a macroscopic volume V mac of the whole 
material, whereas the energy for the formation of a crack 
is considered to be localized in some local volume Vi oc in 
the vicinity of the crack. Let us now imagine the initial- 
ization of an infinitesimally small crack into the solid at 
equilibrium. This perturbation may either open a finite 
crack or the solid may respond in a purely elastic way by 
distributing the energy uniformly over its macroscopic di- 
mensions, which constitutes some kind of unstable equi- 
librium. Then, one can set equal the elastic response 
and the crack formation for infinitesimally small crack 
sizes. Consequently, a correlation between the energies 
(both, elastic as well as cleavage energies) in the localized 
volume Vioc, or -optionally- in the non-local macroscopic 
volume V m ac should exist. This concept of an unstable 



equilibrium should be suitable for general modes of crack 
formation and elastic responses. It might, however, be 
very difficult to derive a correlation in analytic form. 

At least for simple crack modes an analytic formula- 
tions exists. We demonstrate this for cracks of mode I 
as sketched in Fig. ^ along a plane with fixed area A 
and orientation [hkl] a single crystal is cleaved into two 
blocks. The local volume V/ oc can then be expressed by 



Vic 



A L, 



(1) 



with L being a finite length in which the energy is lo- 
calized. Consequently, the macroscopic volume V mac is 
defined as 



Vrr, 



A D, 



(2) 



where D describes the macroscopic thickness of the ma- 
terial. By Eqs. ^ an d0 a simple rescaling condition be- 
tween the two volumes exists, namely V mac = Vi oc (D/L). 
The rescaling factor D/L (or its inverse) transforms lo- 
cal quantities into non-local ones (or vice versa for the 
inverse factor). We will make use of the conversion factor 
when we consider two extreme models of crack formation: 
a) the ideal brittle case modelling cleavage of the mate- 
rial into two rigid blocks, i.e. the atomic positions within 
the blocks do not relax. This process corresponds to a 
crack formation which is so fast that the atomic posi- 
tions do not relax during the process (see panel b of Fig. 
0); b) a process we will call ideal elastic crack forma- 
tion, which is so slow that atoms will always fully relax 
after a crack with a given size x has been openend. If 
x is larger than a critical value then the material should 
abruptly break into two blocks with geometrically relaxed 
surfaces (as sketched in panels c and d of Fig. For 
both cases, analytical models for the crack formation en- 
ergies can be formulated. The yet unknown parameters 
of these models will be determined by fitting the ana- 
lytic expressions to energies obtained from density func- 
tional theory (DFT) calculations modelling the processes 
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a) and b). For this purpose we make use of the Vienna 
Ab initio Simulation Package |!j for calculating the to- 
tal DFT energies and related quantities for suitably large 
slabs. Atomic relaxations are determined by minimizing 
the forces acting on the atoms. 

The final goal, the correlation between elastic proper- 
ties and crack formation, is then formulated for infinitesi- 
mally small crack sizes and elastic strains by setting equal 
crack formation energy and elastic energy, and choosing 
either a local or non-local description. In order to demon- 
strate the generality of our concept we apply it to a vari- 
ety of materials representing different types of bonding, 
such as metallic, covalent and ionic bonding. 



2. IDEAL BRITTLE CLEAVAGE 

As an analytical model for ideal brittle cleavage we 
choose the so-called universal binding relation (UBER) 
|lf)j | for the energy E b (x) depending on the crack size x, 



E b (x) = G b 



x \ fx 

1+ d cxp U 



(3) 



which represents some kind of interatomic potential. The 
material and direction dependent parameters are the 
cleavage energy G b and the critical length l b at which 
the stress o~(x) = dE/dx reaches a maximum (the index 
b denotes brittle cleavage). The UBER energy behaves 
quadratically at x = 0, and for very large separations 
x it approaches the energy zero with a negative curva- 
ture. It should, however, be mentioned that the UBER 
is not universal as usually claimed (see e.g. Refs. 0,E3) 
because the curvature for E b [x — > oo) could also be pos- 
itive [TT| . Furthermore, the two material parameters G b 
and l b might also depend on the actual planes between 
which the crack is initialized if there are different layer- 
to-layer distances along the direction [hkl] (e.g. for the 
diamond lattice in [111] direction). 

In the case of ideal brittle fracture the atomic positions 
remain fixed and do not relax. The electronic states, 
however, will change due to the change of the atomic 
environment in the neighbourhood of the crack because 
potentials at the surface will be different compared to 
the bulk environment. Directional bonds might be re- 
ordered which is reflected in the change of the electron 
distribution and, consequently, the energy of the system. 
All the electronic relaxation effects are fully taken into ac- 
count by the self-consistency of the DFT approach which 
is used to fit the UBER. For this purpose we calculate 
DFT total energies -Ed_ft(^) for a suitably large number 
of different crack sizes Xi and for sufficiently large slabs. 

Studying UBER for a variety of different materials 
with different types of bonding, we found that for most 
cases the fit to the DFT values is rather good and there- 
fore the fitted parameters G b and l b listed in Table |U are 
physically meaningful. Deviations are found for crack 
sizes x > l b which could be particularly large when strong 
covalent bonds are broken. Anyway, the shortcomings of 



UBER are not of significant importance for the present 
study. 

If the response to the perturbation is elastic, the energy 
of the perturbation is expressed by 



E e i ast {x) = ^CVS 2 , 



(4) 



which comprises the uniaxial modulus C for straining 
the material along [hkl] with a fixed area A of the planes 
perpendicular to [hkl] . The symbol V denotes the volume 
over which the energy is uniformcly distributed, and S is 
the dimensionless relative strain. If we decide to refer 
to a local description then V is substituted by V; oc , and 
the relative strain is defined by 5i oc — x/L b with x being 
the initial very small crack size. The localization length 
for brittle cleavage L b which is an intrinsic parameter, 
depends on the material, direction and cleavage planes 
like G b and l b . For very small x we set now equal the 
elastic energy E e i ast and the UBER expression (i.e. the 
T aylo r expansion of E b (x) in Eq. in lowest order of 
x [lj|), which results in the relation 



1 x 2 1 x 2 



(5) 



This is the key equation to correlate the ideal brittle 
cleavage and its material parameters G b ,l bl and L b to 
the elastic properties described by the uniaxial modulus 
C. From Eq. a variety of relations can be derived. For 
the critical stress, i.e. the maximum a b := o~(x = l b ), one 
obtains the equation 

Q~b _ }_^b_}_ _ li^£< (Q\ 

A e A l b e L b 



3. IDEAL ELASTIC CLEAVAGE 

In contrast to the ideal brittle case we now consider a 
cleavage process for which the crack formation is so slow 
that -after opening a crack of size x- the atoms have suf- 
ficient time to fully relax. Furthermore, we require that 
the material should behave perfectly elastic up to a crit- 
ical crack size x = l e , at which it breaks abruptly with- 
out any further interaction between the cleaved blocks, 
as sketched in Fig. n(P ane ls c an d d). We will call this 
process ideal elastic cleavage and we characterize the cor- 
responding parameters by the index e. Denoting the cor- 
responding cleavage energy by G e we derive the obvious 
equation 



E e {x) 



Ge 
11 



(7) 



for the corresponding cleavage energy E e for crack sizes 
x < l e . For crack sizes x > l e we require the condition 
E e (x) = G e . Of course, this is an idealized model because 
a real material will deviate from the simple ideal elastic 
behaviour at least close to the critical crack size, as can 
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be seen in the right panel of Fig. [5]in which the analytic 
model is compared to a "real" material as calculated by 
the DFT procedure. 

Again, like for the brittle case we establish the corre- 
lation between elastic and cleavage properties by setting 
equal elastic and cleavage energy for very small crack 
sizes x. In order to demonstrate the rescaling feature 
of our concept (as discussed at the beginning) we want 
a formulation in terms of a non-local description: the 
strains are now defined with respect to the macroscopic 
thickness by 5 — x/D as in standard elastic theory. Do- 
ing that, we have to rescale the cleavage energy by the 
obvious factor L e /D. Now we have again introduced a 
localisation length as a new material parameter, and we 
can derive the key equation for correlating elastic and 
cleavage properties for the ideal elastic case, 



G - w = l ADC m 



(8) 



It should be noted that in this equation the macroscopic 
thickness D cancels out and, because of this cancellation, 
the relation contains only intrinsic material parameters. 
In fact, we would have obtained the same relation if the 
strain would be defined relative to the localization length 
by 5 — x/L e , similar to treatment of the ideal brittle 
cleavage. 

Exploiting Eq.[S]we derive relations involving the crit- 
ical stress <7 e := a(x = l e ), 



A ~ A L ~ L P 



(9) 



The very close formal relationship to Eq. EJcan clearly 
be seen. However, it should be noted that the material 
dependent cleavage quantities such as G, / and L are dif- 
ferent for both cases (as denoted by the indices b and e) , 
and the constant prefactors are also different. 



4. RESULTS AND DISCUSSION 

In Fig. |2]for the [100] cleavage of NiAl the energies and 
stresses derived from the analytic models are compared to 
DFT results. In general for the brittle case (left panels), 
the agreement is good for smaller x whereas for x > li, 
small but characteristic deviations occur (note also the 
remarks in Ref pl|'). 

The ideal elastic cleavage energy (right top panel) also 
fits quite well to the DFT data. Only close to the criti- 
cal crack size x ss l e deviations from the simple analytic 
model are significant because then atomic relaxations at 
the crack surfaces are becoming important. The differ- 
ences are naturally more pronounced for the stress, be- 
cause it is the first derivative of the cleavage energy. 



Studying the results for the new material parameter, 
the localisation length L, we find a rather striking prop- 
erty for the brittle cleavage as shown in Fig. [3] and in 
Table [J namely that -within an error of ~ 20%- the lo- 
calization lengths, Lb sw 2.4 A , are rather constant for 
all studied materials independent of the directions [hkl] 
and the layer separations oq. Based on this finding we 
derive from Eq. |S|the simple relation 



4.21 



A 



C, 



(10) 



in which cleavage properties (left side) are separated from 
elastic properties (right side) . This equation might serve 
as a starting point for the development of an " engineer- 
ing" tool to estimate the critical cleavage properties in 
terms of the uniaxial modulus C which is a linear com- 
bination of the elastic constants 0, ^1 • 

For the ideal elastic cleavage the critical lengths l e and 
localization lengths L e are much larger than for the brit- 
tle shown in TablellTland Fig-El This seems to be 
obvious because for the ideal elastic cleavage the material 
is now allowed to relax after the crack initialization and 
therefore it needs much larger crack sizes to break it. We 
notice also a strong variation of L e which is in contrast 
to the brittle case. Some -but no simple- correlation be- 
tween l e and L e exists because, generally, for larger l e the 
values of L e are also larger. More studies are needed in 
order to -hopefully- find useful correlations between the 
material parameters for the ideal elastic cleavage. Also 
the critical strengths a e are significantly enhanced from 
20% to 100% in comparions to a^, whereas the cleavage 
energies G e -although reduced compared to G&- differ not 
very strongly from the ideal brittle case (with the excep- 
tion of VC [100]). 

Summarizing, we presented a general concept, which 
correlates elastic and cleavage properties of a solid. By 
its application to two different models for cracks of mode 
I analytic and quantitative correlations were derived by 
introducing a localization length L , a new intrinsic mate- 
rial parameter which proved to be meaningful for a large 
variety of materials. Future efforts should focus on the 
fundamental understanding of L. 
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Tables 



TABLE I: Parameters for brittle cleavage for some selected 
compounds and directions: uniaxial elastic modulus C (GPa), 
cleavage energy per surface area Gb/A (J/m 2 ), critical length 
lb (A ), maximum stress o~b/A (GPa), interlayer distance be- 
tween cleaved planes oo (A), and localization length Lb (A 
) for selected compounds and cleavage directions [hkl]. In 
brackets, the crystal structure is denoted in terms of Pear- 
son's notation. 







[hkl] 


c 


Gb/A 


lb 


o-b/A 


ao 


L b 


NiAl 


(B2) 


100 


203 


4.8 


0.69 


26 


1.45 


2.0 






110 


284 


3.2 


0.54 


22 


2.05 


2.5 






111 


311 


4.1 


0.58 


26 


0.84 


2.4 


Ni 3 Al 


(Ll a ) 


100 


225 


4.3 


0.66 


24 


1.78 


2.1 






111 


331 


3.7 


0.52 


26 


2.06 


2.4 


VC 


(Bl) 


100 


647 


3.2 


0.37 


32 


2.08 


2.8 






110 


585 


7.0 


0.55 


46 


1.47 


2.5 






111 


564 


9.9 


0.58 


63 


1.20 


1.9 


MgO 


(Bl) 


100 


299 


1.8 


0.37 


18 


2.11 


2.2 






110 


345 


4.4 


0.54 


30 


1.22 


2.3 



TABLE II: 
Table □ 



Ideal elastic cleavage: Similar quantities as in 





[hkl] 


c 


G e /A 


le 


O-e/A 


L e 


NiAl 


100 


203 


4.6 


2.7 


34 


15.8 




110 


284 


3.1 


2.0 


32 


17.7 




111 


311 


3.9 


2.2 


36 


18.4 


Ni 3 Al 


100 


225 


3.0 


2.2 


27 


18.2 




111 


331 


3.2 


1.0 


64 


5.2 


VC 


100 


647 


2.4 


0.8 


60 


6.5 




110 


585 


6.0 


1.6 


75 


12.5 




111 


564 


8.4 


1.6 


105 


8.6 


MgO 


100 


299 


1.7 


0.8 


42 


5.3 




111 


345 


10.2 


1.9 


107 


6.9 
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FIG. 1: Models for cracks of mode I. A solid (sketched as a 
stacking of interacting layers with a layer distance ao, panel 
a) undergoes brittle cleavage as sketched in panel b: a crack 
of size x cleaves the material into two rigid blocks without 
relaxing the atomic geometry of the layers in the blocks. For 
the ideal elastic cleavage (see text) the material relaxes in a 
perfectly elastic manner (panel c) up to a critical crack size 
above which it breaks abruptly into two blocks with atomi- 
cally relaxed surfaces (panel d). 
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FIG. 2: Cleavage along [100] of NiAl: energy and stress 
per surface area vs. crack size x for the ideal brittle (left 
panel) and ideal elastic (right panel) cleavage model full lines: 
analytic models, dots: DFT results). 



FIG. 3: Localization lengths vs. critical lengths for ideal 
brittle and elastic cleavage for a variety of materials and di- 
rections. 




